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' Phantom fields with negative kinetic energy are often plagued by the vacuum quantum instability 

in the ultraviolet region. We present a Lorentz-violating dark energy model free from this problem 
and show that the crossing of the cosmological constant boundary w = — 1 to the phantom equation 
of state is realized before reaching a de Sitter attractor. Another interesting feature is a peculiar time- 
dependence of the effective Newton's constant; the magnitude of this effect is naturally small but may 
be close to experimental limits. We also derive momentum scales of instabilities at which tachyons 
or ghosts appear in the infrared region around the present Hubble scale and clarify the conditions 
under which tachyonic instabilities do not spoil homogeneity of the present/future Universe. 
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The compilations of various observational data show that the Universe has entered the stage of an accelerated 
expansion around the redshift z ~ 1 [H 0, H, 0, S @] • The equation of state (EOS) parameter w of Dark Energy (DE) 
responsible for the acceleration of the Universe has been constrained to be close to w — — 1. However, the phantom 
EOS (w < —1) is still allowed by observations and even favored by some analyses of the data Q. It is also possible 
that the EOS of DE crossed the cosmological constant boundary (w = — 1) in relatively near past Q. 

The presence of the phantom corresponds to the violation of weak energy condition, the property which is generally 
difficult to accommodate within the framework of field theory. The simplest model which realizes the phantom EOS is 
provided by a minimally coupled scalar field with a negative kinetic term [!G3 (see also Refs. [HlGj). The negative 
kinetic energy is generally problematic because it leads to a quantum instability of the vacuum in the ultraviolet 
t^. | (UV) region 

M HE Hi HE 111 HI: the vacuum is unstable against the catastrophic particle production of ghosts and 
■ normal (positive energy) fields. 

There have been a number of attempts to realize the phantom EOS without having the pathological behaviour in 
• . the UV region. One example is scalar-tensor gravity in which a scalar field <f> with a positive kinetic term is coupled to 
Ricci scalar R This coupling leads to the modification of gravitational constant, but it was shown in Ref. [2(| 

that there are some parameter regions in which a phantom effective EOS is achieved without violating local gravity 
constraints in the present Universe. 

Another example isprovided by the so-called modified gravity, including f(R) gravity models [2lJ and the Gauss- 
Bonnet (GB) models [22|]. In f(R) models it is possible to obtain a strongly phantom effective EOS, but in that case 
the preceding matter epoch is practically absent [23| . For GB DE models it was shown in Ref. [HI that the crossing of 
the cosmological constant boundary, w = —1, is possible, but local gravity experiments place rather strong constraints 
on the effective GB energy fraction (25[. In addition, tensor perturbations are typically plagued by instabilities in 
the UV region if the GB term is responsible for the accelerated expansion of the Universe [2y|. Thus, it is generally 
not so easy to construct viable modified gravity models that realize the phantom effective EOS without violating 
cosmological and local gravity constraints. 

The third example is the Dvali-Gabadadze-Porrati (DGP) braneworld model 27] and its extension [2Sfl with a GB 
term in the bulk, which allow for the possibility to have w < — 1 [1^, [3(|. However, it was shown in Ref. |3l| that the 
DGP model contains a ghost mode, which casts doubts on the viability of the self- accelerating solution. 

While the above models more or less correspond to the modification of gravity, it was recently shown that in the 
Einstein gravity in a Lorentz-violating background the phantom EOS can be achieved without any inconsistency in 
the UV region |3^.[33|. In particular, in the model of Ref. [33| Lorentz invariance is broken in the presence of a vector 
field which has two-derivative kinetic terms similar to those given in Ref. [H| . The effect of the Lorentz violation 
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is quantified by a parameter 5 = B^B 11 /A4 2 , where M. is an UV cut-off scale. In analogy to Ref. [3|| the vector field 
also has one-derivative coupling ed^B^ with a scalar field <&, where e is a small parameter that characterizes an IR 
scale. In the UV region, where the spatial momentum p is much larger than e, ghosts, tachyons and super-luminal 
modes are not present. Meanwhile tachyons or ghosts can appear in the IR region p < e. This is not problematic 
provided that e is close to the present Hubble scale. 

In this paper we apply this Lorentz-violating model to dark energy and study the cosmological dynamics in detail 
in the presence of mass terms in the potential, V = ^m 2 ® 2 — ^M 2 X 2 (where X 2 = B^B^). We show that the 
model has a de Sitter attractor responsible for the late-time acceleration. At early times DE naturally has normal 
EOS with w > — 1, while the phantom EOS can be realized between the matter-dominated era and the final de Sitter 
epoch. We clarify the conditions under which the cosmological constant boundary crossing to the phantom region 
occurs. Interestingly, in a range of parameters this crossing takes place at the epoch when Q m ~ 51de thus making 
the crossing potentially observable. 

Another interesting feature of our model is the time-dependence of the effective Newton's constant. It is naturally 
weak, but may well be comparable with current experimental limits. Moreover, the effective Newton's constant G*(t) 
has a peculiar behaviour correlated with the deviation of w from — 1 . 

We also derive momentum scales of instabilities of perturbations, first in Minkowski spacetime. This is the extension 
of the work [33[ that mainly focused on the case of massless scalar (m = 0) . We show that in the UV region (p e) 
the model does not have any unhealthy states such as ghosts, tachyons or super-luminal modes. In the IR region 
(p < e) tachyons or ghosts appear, depending on the momentum. Finally, we study the evolution of perturbations 
in the cosmological background and estimate the amplitude of perturbations amplified by the tachyonic instability 
around the scale of the present Hubble radius. The perturbations remain to be smaller than the background fields 
under certain restriction on the model parameters. 

This paper is organized as follows. In Sec.|H]we present our Lorentz-violating model and derive basic equations 
describing spatially flat Friedmann-Robertson- Walker cosmology in the presence of DE, radiation and non-relativistic 
matter. In Sec. IIHI the cosmological dynamics is discussed in detail analytically and numerically with an emphasis 
on the occurrence of a phantom phase before reaching a de Sitter attractor. The time-dependence of the effective 
gravitational constant is also considered. In Sec. IIVI we study the Minkowski spectrum of field perturbations and 
clarify the properties of tachyons and ghosts in the IR region. We then discuss the tachyonic amplification of field 
perturbations around the present Hubble scale in the cosmological background. We summarize our results in Sec. [V] 
Appendix A contains the derivation of the effective "Newtonian gravitational constant" in our model. In Appendix 
B we derive the fixed points of the system by rewriting the equations in autonomous form. We analyse the stability 
of the fixed points and show analytically that the cosmological evolution proceeds from radiation-dominated stage 
through matter-dominated stage to the final de Sitter regime. 



II. LORENTZ-VIOLATING MODEL 



We study a 4-dimensional Lorentz-violating model whose Lagrangian density includes a vector field B^ and a scalar 
field $: 

C = ~a(E)g vX D l ,B v D^B x + ±p(Z)D^B„D»B x ^^- + ±d M *d»$ + ed^B" - V(B, *) , (1) 

where S = B^B^ / A4 2 with M. being an UV cut-off scale of the effective theory. The dimensionless parameters a 
and (3 are the functions of S, and e is a free positive parameter that characterizes an IR scale. The first two terms 
in |lj are familiar in two-derivative theory [34|], whereas the one-derivative term ed^B^ is introduced following the 
approach of Ref. [35[ . 

We study dynamics of flat Friedmann-Robertson- Walker (FRW) Universe 

ds 2 = Af 2 {t)dt 2 - a 2 (i)dx 2 , (2) 

where Af(t) is a Lapse function and a(t) is a scale factor. In the case of spatially homogeneous fields with B. L = 
(i = 1,2,3), the Lagrangian |T]) reads 

V~gC = Y^X 2 - ^X 2 + ~^ 2 + ea^X - a*MV(X, 0) , (3) 
where X = Bq/M, <j> is the homogeneous part of the field $ and 



-,iX) = ^0(X)-a{X). (4) 
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Hereafter we study the case in which the following condition holds 

a > 7 > . 

This is required to avoid a superluminal propagation in Minkowski spacetime [33| . as we will see later. Throughout 
this paper we assume that a and 7 are of order unity. 

For fixed X, the second term in the Lagrangian ([3]) has precisely the form of the Einstein-Hilbert action specified 
to the flat FRW metric. Hence, it leads to the change of the "cosmological" effective Planck mass |33[ 

mpl,cosm = TO pi + 47raX 2 . (5) 

Another effective Planck mass ra p i,Newton determines the strength of gravitational interactions at distances much 
shorter than the cosmological scale; in general, these two effective Planck masses are different fl9l. [35l l36l|. We show 
in Appendix A that the "Newtonian" Planck mass in our model is given by 

m P LNcwton = m li - 47raA 2 . (6) 

Both effective Planck masses depend on time via X — X(t). Since the time-dependent terms in ([5]) and ([6]) differ by 
sign only, it will be sufficient to study one of these effective masses. In what follows we concentrate on the "Newtonian" 
mass © for definitencss. 

In this paper we focus on the case in which the potential V takes a separable form: 

V = W{<t>) + U(X) . (7) 

We take into account the contributions of non-relativistic matter and radiation whose energy densities p m and p r , 
respectively, satisfy 

pm + 3Hp m = , (8) 

pr + AHp r = . (9) 

The energy density of the fields is derived by taking the derivative with respect to Af of the action S = J d 4 Xy/—g£: 

SS 



I 

P = 3 



SAf 



= - ^-H 2 X 2 + \tf + V. (10) 



We set Af = 1 for the rest of this paper. 
The Friedmann equation is given by 



H 2 = 



, 2 

a 



\px 2 - ^h 2 x 2 + x -tf + W{4>) + U(X) + p m + Pr 



(11) 



where k 2 = 877/771^. The equations of motion for the homogeneous fields <f> and x are 

--y(jt + 3#X) - \l,xX 2 - ^a, x H 2 X 2 - 3aH 2 X + e<j) = U x , (12) 

—((f) + 3H<fi) - e(X + 3HX) = W A , (13) 
where 7.x = d7/dA, etc. Taking the time-derivative of Eq. (jlip and using Eqs. (TT2|) and p^| . we obtain 

ij k2 f 4 

H = ~— Ip+P + Pm + -Pi 

where 

p + p = e<t>X + aHX 2 + 2aHXX + 7A 2 + <j) 2 + a, x HX 2 X . (14) 

In what follows we assume for simplicity that a and 7 are constants, i.e., a t x — 7.x = 0. 
Following Ref. [H| we consider the simplest potential for the fields, 

W(<t>) = \m 2 4> 2 , U(X) = - l -M 2 X 2 , (15) 
which allows for a possibility to realize a phantom phase. 
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III. DYNAMICS OF DARK ENERGY 



One way to analyse the cosmological dynamics in our model is to make use of the autonomous equations, the 
techniques widely used in the context of dark energy studies [1, [H, This approach is presented in Appendix B, 
where we analytically confirm that our model can lead to the sequence of radiation, matter and accelerated epochs. 
Also, in Appendix B we derive the conditions under which the de Sitter solution given below is an attractor. Here we 
first present a simpler analysis based on the slow- roll approximation. Then we give numerical solutions to eqs. ([H]), 
©, (|lip . (|12jl . (|13p . exhibiting transient phantom behaviour, and study their dependence on various parameters of 
our model, including the initial values of the fields. 



A. Final and initial stages 

One immediate point to note is that in the absence of radiation and matter, the system of equations (112p . 
T3| has a de Sitter solution, H = const, for which <f> and X are also independent of time, provided that 

e 1 2a , „ 

->\hr- 16 

m V o 



Indeed, for constant H, <fi and X eqs. (fTTj) , (TT2"]) . (| 13|) reduce to a simple algebraic system 

H 2 
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.^H*X>-^X> + ^ 
2 2 2 v 



3aH 2 = M 2 , 

-3eHX = m 2 cj). (17) 
Once the inequality §TE§ is satisfied, this system has a solution 

Ha - M 



3a 

3 Mm p ie 



4tt y/a~m 2 ^fWjm? - 2a ' 

A A = -^ 1 (18) 

V47T \Jie 2 /m 2 — 2a 

We will see in what follows, and elaborate in Appendix B, that in a range of parameters this solution is an attractor 
which corresponds to the de Sitter phase in asymptotic future (hence the notation). In order to use this for dark 
energy we require that the mass scale M is of the order of the present Hubble parameter H . Then the Newtonian 
effective Planck mass, Eq. (j6]), is given by 

™pl,Ncwton = Wpl (l " 3e 2 /m a 2_2a) ' (19) 

In order that the change of the Planck mass be small, we impose the condition 

e > Vara . (20) 

It is worth noting that under this condition, the contribution of the field <p in the energy density dominates in the 
de Sitter regime, 

m 2 ,, M 2 „ 3aH 2 „ . . 

—4> 2 A » — X\ - —j-xl . (21) 

Thus, as the system approaches the de Sitter attractor, the total energy density in the Universe becomes determined 
by the scalar field energy density. 

Another point to note is that at early times (at the radiation-dominated epoch already), when the Hubble parameter 
is large enough, the term (— 3aH 2 X) in Eq. (fT2")) drives the field X to zero, the relevant time being of the order of 
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the Hubble time. Soon after that the field <p obeys the usual scalar field equation in the expanding Universe, so the 
Hubble friction freezes this field out. Thus, the initial data for the interesting part of the DE evolution are 

4>i = const . (22) 

The value of <j>i is a free parameter of the cosmological evolution in our model. Since at early times the field X is 
close to zero, its effect on the evolution of the field 4> is negligible. The field <f> slowly rolls down its potential, and its 
energy density dominates over that of X. Therefore, EOS for DE at early times is normal, w > — 1, with w being 
close to —1. We refer to this regime as quintessence stage. As we will see below, in a range of parameters, the system 
eventually crosses the cosmological constant boundary w — — 1 and passes through a transient phantom phase before 
reaching the de Sitter asymptotics (p~8|) . 

B. Slow roll phantom regime 

The approach to the de Sitter solution (|18p occurs in the slow roll regime. To see how this happens, we truncate 
Eqs. (TTJ]) and (13]) to 

e<j) - 3aH 2 X = U }X , (23) 
-ZeHX = . (24) 

This truncation is legitimate provided that in addition to the usual slow-roll conditions 4> -€i Hcf> and X <§; HX, the 
following conditions are satisfied: 

0<eA, (25) 
e<t>X <C V , (26) 
X <C HX . (27) 

[When writing inequalities, we always mean the absolute values of the quantities.] Note that we do not impose the 
condition e0 ^> 3aH 2 X unlike in Ref. since the term 3aH 2 X is not necessarily negligible relative to the term 
U.x in Eq. 

From Eq. (f!?3"|) we obtain 



where 



{ = l-3a^. (29) 



Note that £ may be considered as a measure of the deviation from the de Sitter regime (fT5|) . 
Substituting Eq. (|28|) into Eq. ([24]) we get the following equation 



ZH<P - SiW,<p , (30) 

where 

Equation (|30|) shows that the field 4> rolls up the potential W(4>) for £ > 0, i.e., for 

H<^L. (3D 
V3a 

This is the region in which the phantom equation of state (w < — 1) is realized; indeed, Eq. (|14p gives p + p ~ e<j>X — 
£XU x = — £M 2 X 2 . Another way to understand the phantom behaviour is to notice that when the system approaches 
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the de Sitter regime, the field <j> dominates the energy density, see Eq. (f2Tj) . so the energy density increases as the 
field <fi rolls up. 

Let us find out whether the slow roll conditions (|25|) . (I26| and (|27|) are indeed satisfied. Making use of Eq. (|28|) we 
obtain that the condition (f2"5| is equivalent to 

e 2 > £Af 2 , (32) 
while using Eqs. ([28| and (|30|) we rewrite the condition (|26|) as 

2 „m 2 M 2 

e » ^~h5- • ( 33 ) 



The second inequality ensures also the validity of the relation (|2T|) : this can be seen by taking the time derivative 
of Eq. The latter two inequalities are automatically valid at small £, that is near the de Sitter solution (|18p . 

We conclude that the approach to the de Sitter solution indeed occurs in the slow roll regime, and that the phantom 
phase is indeed realised provided that the relation f3~l]) holds. Our analysis implies also that the de Sitter solution 
(I18p is an attractor: the Hubble parameter slowly increases towards its de Sitter value, £ decreases, and the dynamics 
gets frozen as £ — > 0. 

Since the field </> dominates the energy density at the phantom slow roll stage, the condition (f3"Tj) takes a simple 
form 

*<*A = -^_, (34) 

where we made use of (f2T)]) . The latter relation translates into the range of initial conditions which eventually lead to 
the transient phantom behaviour, 



< 



(35) 



Indeed, during the radiation- and matter-dominated stages the field (f> remains almost constant, and at the quintessence 
stage it also does not roll down much. 

Recalling again that <f> dominates the energy density, we rewrite the inequality (I33[) as e s -C 1, where 

2am 2 ( (j) A \ 2 

The parameter e s may be viewed as the slow roll parameter for the field <p. Indeed, one observes that 

Xr- = £ s , 37 

which, together with Eq. (|30|) , justifies this interpretation. 

From Eqs. fTD]). (TH)) . (j2"5|) and $5U§ it follows that during the slow roll phantom stage, the EOM parameter of DE 
is given by 

w = -l-e s . (38) 

Hence, the appreciable deviation from w = — 1 occurs when e s is not much smaller than unity, i.e., when (pi is 
appreciably smaller than tp^. 

In the next section we confirm these expectations by numerical analysis, and also show explicitly that in a range of 
parameters, the cosmological evolution proceeds from radiation-dominated to matter-dominated epoch, and then to 
the slow roll phantom stage, before finally ending up in the de Sitter regime fT8]) . 



C. Numerical solutions 

In our numerical analysis we choose initial conditions X = <f> = X = with nonzero values of <f>, p m and p r . This 
choice corresponds to the initial data (|2"2"|) . We have also tried many other initial conditions and found that the results 
are not sensitive to the initial values of X, X and </>, in accord with the discussion in the end of Sec. IIII Al 
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In Fig. [T] we plot the cosmological evolution for the model parameters a = 1, 7 = 1/2, e/m 
initial value <f>i — 0.5(f>A- We find it convenient to present the plots in terms of the variable 



3, M/m = 1 and the 



N = baa. 



(39) 



Figure [T] clearly shows that the sequence of radiation, matter and de Sitter epochs can be achieved together with the 
w = — 1 crossing. The DE EOS parameter w is nearly a constant, w ~ — 1, during the radiation and matter epochs 
because the fields are almost frozen. At the transition era from matter domination to DE domination, w begins 
to grow because the kinetic energies of the fields become important; this is the quintessence phase. However, the 
system soon enters the phantom phase during which the field <$> rolls up the potential. Hence the equation of state w 
crosses the cosmological constant boundary w = — 1 and reaches a minimum value u> m in < — 1- The solution finally 
approaches the de Sitter regime (|18[) from the phantom side. Of some interest is also the EOS parameter of the entire 
system, 



Wcff 



Ptot 
Ptot 



(40) 



It is seen from Fig. [TJ this parameter also becomes smaller than —1 soon after w = — 1 crossing. 
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FIG. 1: Cosmological evolution for the model parameters a = 1, 7 = 1/2, e/m = 3 and M/m = 1. We choose initial conditions 
X = X = (j> ~ 0, </>i = 0-5(f>A, and f2 r ,i = 0.99, Q m ,i = 0.01. Shown is the evolution of Ode, f2m, Sir, w and u> e fr as functions 
of N = In a. Note that the present epoch corresponds to Q. m — 0.3 and Ode = 0.7, which is denoted by a vertical line. After 
the cosmological constant boundary crossing, the DE EOS parameter w reaches a minimum w m in = —1-19 and then increases 
towards the de Sitter value w = —1 from the phantom side. 



The pattern shown in Fig. |T| is generic in our model, provided that its parameters and initial data obey m/e <C 1, 
M/e <§; 1 and <j>i < 4>a (in fact, the inequalities here need not be strong). The strengths of the effects depend, of 
course, on the values of these parameters. In particular, the minimum value u> m i n is related to the slow-roll parameter 
e s , in accord with Eq. (|38|) . If the field cf> evolves very slowly, one has e s -C 1, so uWn is close to —1. On the contrary, 
the appreciable deviation from w = — 1 occurs if e s is not very much smaller than unity. 
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FIG. 2: The minimum value of the EOS parameter w of DE, its value wo at the present epoch (fl m = 0.3, Qbe = 0.7) and its 
maximum value as functions of m/e for M/e = 1/30 and the initial value <j>i = 0.50a- 



Since the field <fi is practically frozen during the radiation- and matter-dominated epochs, and evolves rather slowly 
later on, the dependence of the cosmological evolution on the parameters of the model and on the initial value <fii can 
be understood, at qualitative level, by inspecting Eq. ([55)) . For the qualitative discussion of the evolution well before 
the asymptotic de Sitter regime sets in, the parameter £ in (f36|) may be set equal to 1, while the value of <fi may be set 
equal to its initial value <f>i. Then Eq. (|36p implies that with other parameters and <f>i fixed, for smaller m/e one gets 
smaller absolute value of e s at its minimum, leading to the value of u> m i n closer to —1. This is shown in Fig. [3J From 
Fig. [1] it is clear, however, that the minimum of w occurs after the present epoch (£!de — 0.7 and Q m ~ 0.3). Again, 
this is a rather generic feature of our model. Therefore, instead of w m in, more interesting quantities are the present 
value wo of the DE EOS parameter and also its maximum value before the cosmological constant boundary crossing. 
These quantities are also given in Fig. [5] Overall, the behaviour shown in Fig. [T]is more pronounced at larger m/e, 
once other parameters of solutions are kept fixed. 

According to Eq. (|3"f5|) , the initial value of the field 4> is also important to determine the amplitude of e s and hence 
w: a smaller 4>i results in a stronger deviation of w m { n from —1. The present value wq also becomes more negative 
(stronger deviating from —1), while the cosmological constant boundary crossing occurs earlier. From the numerical 
analysis we find that at relatively large values of m/e the increase of M has the opposite effect, while at smaller m/e 
the effects due to the variation of M are small. These properties are illustrated in Figs. [3] and [U 

As discussed above, the phantom phase occurs only if the initial value of <j> obeys (|35p . otherwise the de Sitter 
attractor is approached from the quintessence side, w > —1. This is illustrated in Fig. [5] 
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FIG. 3: The dependence of the present value wo (i.e., w at Q m = 0.3, Ode = 0.7) on the initial value of <f>, for different sets of 
the model parameters. 
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FIG. 4: The dependence of the redshift (from the present epoch, Q. m — 0.3, Ode = 0.7) of the cosmological constant boundary 
crossing, w — —1, on the initial value of (f>. 
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FIG. 5: Cosmological evolution for the model parameters a = 1, 7 = 1/2, e/ra = 3 and M/m = 0.1. We choose initial 
conditions X — X = <f> — 0, 4>i — 3(j>A, and £l r< i — 0.99, fl m ,i = 0.01. The present epoch corresponds to Qde = 0.7 and 
fl m = 0.3, which is denoted by a vertical line. In this case the cosmological constant boundary crossing is not realized because 
<fi remains always larger than (f>A- 

At the end of this section we discuss the variation of the effective gravitational constant. According to Eq. |6]), the 
effective Newton's constant that determines the interaction deep inside the horizon scale is given by 



G* = G(l - inGaX 2 )- 1 . 

Its variation in time is conveniently expressed in terms of the following quantity 

dlnG* _ G* _ 8naXX 
din a = G~H ~ H{m 2 pl - 4naX 2 ) ' 



(41) 



(42) 



The typical experimental and observational constraints on the variation of G* in the present Universe are given by 
|G*/G*| < 10~ 12 yr -1 [37j], which translates into the condition 



G*/G* < 10" 2 i? . 



(43) 
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FIG. 6: The evolution of the variation of the effective gravitational constant d In G. /din a as function of N = In a and redshift 
z. The cases (a) and (b) correspond to the model parameters and initial conditions given of Figs. [T]and[5] respectively. The 
black points represent the values at the present epoch (fi m ~ 0.3). 
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FIG. 7: The maximum value of din G, /din a as a function of m/e. Other parameters and initial conditions are the same as in 
Fig. El 



In Fig. [5] we plot the evolution of the quantity d In G* /din a for the model parameters and initial conditions given in 
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Figs. [T]and [5] At the present epoch (Vt m ~ 0.3) we obtain the values G*/G* = 3.5 x lO^H"^ 1 and 2.4 x 10~ 2 H ( J 1 for 
these two cases, respectively. Comparing Fig. [5] with Figs. [T] and [5] one observes that the variation of the gravitational 
constant is correlated in time with the deviation of w from — 1. This is clear from Eq. (|4"Tj) too: the gravitational 
constant varies when the field X changes in time, while the latter occurs during the transition from the matter 
dominated stage to the final de Sitter attractor. It is precisely at this transition stage that w substantially deviates 
from — 1. 

Figure [7] shows the maximum value of d In G* /din a as a function of m/e. Again, the variation of the effective 
Newton's constant is more pronounced at larger m/e. This means that it is correlated with the amplitude of the 
deviation of w from —1. The dependence of d In G* /din a on M and on the initial value of (j> is rather weak. 

It is worth pointing out that as long as the deviation of the EOS from w = — 1 is not so significant, the models 
satisfy the constraint (|4*3"|) , and also that our model suggests that the variation of G» is close to the present upper 
bound on G»/G*. 



IV. MOMENTUM SCALES OF INSTABILITIES 



In this section the momentum scales of instabilities are present in our model. We first study dispersion relations in 
the Minkowski space-time and then proceed to those in the FRW space-time. We wish to clarify the conditions under 
which a tachyon or a superluminal mode appears by considering dispersion relations. We also evaluate the energy of 
the modes to find out a ghost state. 



A. Minkowski spectrum 



Let us consider the perturbations for the fields, 

B = X + b , Bi = h, $ = cj) + tp . (44) 

The quadratic Lagrangian for perturbations, following from the general expression ([TJ, is 

y ol 1 

£b M,v = -d^bod^bo + -df+bid^bi + -d^ipd^ip + ed (pb - ediipbi 



2 



-mlbl - -m\b\ - -m%<p 2 , (45) 



where 



ml = u ,xx, m i = _ ~x^' m l^ W ,4"t>- ( 46 ) 

For our model ([15)) one has —m\ = m\ = M 2 and m 2 = to 2 . In what follows we concentrate on this case, and assume 
the following relations, see (|2T)|) and (|3"2"1) . 

e > \fam , 

e > M . (47) 

Varying the Lagrangian (|45|) with respect to bi, bg and ip, we obtain the equations for the field perturbations. In 
order to find the spectrum of the system we write the solutions in the form bo = &oe ^p^ ' x, = boe^ UJt ~ p ' r \ bi = bie %v ^ 
and if = ipe tp ^ x . The transverse mode of the vector field Bi has the dispersion relation 

-oW + — ■ (48) 
a 

The three scalar modes bi — {j>i/p)bi,, bo and tp satisfy the following equations 

M 2 \~ e _ „ 

















-P 2 + 




-P 2 -' 



(49) 

a J a 
M 2 \~ e 

&o-?-w^ = 0. 50) 

7 / 7 

to 2 ) tp — ieubo — iepbh = . (51) 
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Expressing 6^ and bo in terms of (p from Eqs. (|49f . (I50|) and plugging them into Eq. (|5T|) . we find that the eigenfre- 
qucncies corresponding to three mixed states satisfy 

Z ' p2 ' p2 ^=0, (52) 





M 2 \ 










(-£) 






'—J 







where 



77a 7a 



a/ - p 2 . (53) 



The spectrum in the case m = was studied in Ref. [33| . Our purpose here is to extend the analysis to the case of 
non-zero m. Denoting the solutions of Eq. (|52l as z±, z 2 and Z3, we obtain the relation 

m 2 M 4 

ziz 2 z 3 = . (54) 

Once the conditions (j47|) are satisfied, then one can show that if the relation 

Z\ < z 2 < z 3 (55) 

holds at some momentum, then the inequality (|55p is satisfied for all momenta. 

In the limits p — > 00 and p — > 0, we obtain the following dispersion relations, respectively. 

• (A) UV limit (p ->■ 00) 



e /l 1 e 2 /2rn 2 11.,., 

= p- + - + — + +O M J /p), 56 

2 V 7 a 8p \ 7 a/ 

^2 = P+ 3 2( - > +Q 1/p 5 ), 
2p' :i e^(a + 7) 



e /I F e 2 /2to 2 1 1\ , , s 

^3 = P+~\ - + - + -£z[— + )+0(M 2 p). 

2 V 7 a 8p \ e z 7 a J 

We see that lo\ < uj 2 < 0J3 and z\ < 0, z 2 ,3 > 0. In all three cases the group velocities duji/dp are less than 1, 
so neither mode is superluminal at high three-momenta, provided that a > 7. The two-derivative terms in the 
Lagrangian (|4"S"1) dominate in the UV limit, so there are neither ghosts nor tachyons in this limit. 

(B) IR limit (p -> 0) 

Wl = -^^[l + 0(m 2 / £ 2 ,M 2 / e 2 )] , (57) 

Af 2 

w 2 = , 

a 

u> 3 = -+m 2 + 0(M 2 ). 

7 

We see again that io\ < lu 2 < ^3 and z\ < 0, 2:2,3 > 0. This means that using the property (|55[) we can identify 
the modes: the first one has the behaviour (l56l) and (l57l). and so on. 



It follows from Eq. (f54|) that z< never vanish. In fact the coefficients in Eq. (|52|) are regular at all momenta, so Zi 
are regular as well. Therefore, z% never change signs and hence z± < 0, Z2.3 > for all momenta. This means, in 
particular, that the second and third modes never become tachyonic. 

Let us discuss the dangerous mode with the dispersion relation uj = lo\(p) in some detail. The expression for the 
fields in each mode is 

h,i = -iep {izi + M 2 ) ■ d , 
b ,i = ieuj(az l - M 2 ) • C 2 , 

ipi = ( 1 z l + M 2 )(az l - M 2 )-d, (58) 
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where Ci are the normalization factors. Setting lo 2 — in Eq. ([52]) . we obtain three critical momenta 

e 2 - M 2 



Pl,2 



i 

~ 2 

M 2 
7 



^ 2 ± 



e 2 -M 2 



2 _ 



4to 2 M 2 



n 



Under the conditions (|47p , the critical momenta p 2 2 are approximately given by 

- 2 M 2 „ n m 2 M 2 



Pi 



— m 



Pi 



so that p\ > p\ > p 2 > 0. The tachyonic mode (ui 2 < 0) is present for < p 2 < p\ and p\ < p 2 < p\ . 
In order to find whether there are ghosts we calculate the energy of the modes ^ v ' 



(59) 
(60) 

(61) 



Eiip) = 2oj 2 \C,\ 2 [aeVfrzi + M 2 ) 2 + e 2 ( 7 p 2 ~ M 2 )(a Zl - M 2 ) 2 + ( 7 z, + M 2 ) 2 {az t - M 2 ) 2 ] . (62) 

For the modes with W2,3 we have -Z?2,3(p) > 0. For the mode with uj\ the energy is equal to zero at p — p% = M/^Fy. 
While uj\ > for p 2 < p\ = M 2 /j, the energy E\{p) changes its sign at this momentum. Thus the mode with u>\ is a 
ghost for p 2 < M 2 /^f. 

We summarize the properties of the dangerous mode as follows: 

• (i) p 2 > (e 2 - M 2 )/a - m 2 : healthy 

• (ii) Af 2 /7 < p 2 < (e 2 - M 2 )/a - m 2 : tachyon 

• (iii) m, 2 M 2 /e 2 < p 2 < M 2 /^: ghost, but not tachyon 

• (iv) < p 2 < m 2 M 2 /e 2 : tachyon. 

Unlike the case m = [33| the tachyon is present in the deep IR region (iv). 

To end up the discussion of the modes in Minkowski space-time, we give the expressions for the minimum values of 



lu 2 in the tachyonic regions 



7e 2 2 e 2 7 + 2a 

~~a — 7 \ at p = 

4a(a + 7 ) 4a 7 + a 



(63) 



(iv): 



2 

w min — 



% 2 M 2 



at p 2 = - 



Note that |wjjj n | is relatively large in the region (ii), so this region is the most problematic. 



B. Evolution of perturbations in cosmological background 

Finally we discuss the evolution of field perturbations in the FRW background In the cosmological context 
the physical momentum p is related to the comoving momentum k as p = k/a. Once the parameters of the model 
and initial data are such that the cosmological boundary crossing occurs, the present epoch (JIde — 0.7) typically 
corresponds to the phantom region. From Fig. 1 one can see that the Hubble parameter does not change much during 
the transition from the phantom epoch to the final de Sitter era. Hence the present value of the Hubble parameter (Hq) 
is of the same order as the value H = M/y/3a in the de Sitter asymptotics. This means that the value P3 = MJ y/j 
is of the same order as Hq provided that 7 and a are of order unity. 

The tachyon appears when the momentum p = k/a of the dangerous mode becomes smaller than 
\J (e 2 — M 2 )/ a — m 2 and temporally disappears when the mode crosses the value M / 4/7. Hence this instability 
is present for the modes which are inside the Hubble radius and satisfy M/^/7 < p < \J (e 2 — M 2 )/a — m 2 , but it 
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is absent for the modes deep inside the Hubble radius, satisfying p > \J (e 2 — M 2 )/a — m 2 . After the Hubble radius 
crossing (k = aH), the tachyonic instability disappears in the momentum region m 2 M 2 /e 2 < p 2 < M 2 /j, but the 
tachyon appears again for p 2 < m 2 M 2 /e 2 . Note that the ghost existing at m 2 M 2 /e 2 < p 2 < M 2 /j is a not a problem 
because of its low energy [l3|, HO] . 

In what follows we discuss the evolution of field perturbations in the two tachyonic regimes. Before doing that it is 
instructive to study the high-momentum regime that sets the initial data for the tachyonic evolution. 

1. p 2 > (e 2 - M 2 )/a~m 2 

We denote the overall amplitude of the dangerous mode as <p. Since the modes are deep inside the Hubble radius 
(k/a 3> H) in the regime we discuss here, the field perturbation x approximately satisfies 



i 2 

d?/ 



■j-zX + k'x-O- 



(64) 



where x = al P an d f] is conformal time defined by r\ = J a dt. Taking the asymptotic Minkowski vacuum state, 
X = e~ lkn /y/2k, the squared amplitude of the field perturbation ip is given by [4l[ 

2 



Ink 6 



m 



k 
2ira 



(65) 



Since the maximum momentum at which the tachyon appears is k/a ~ e/y/a, one has the following estimate for the 
amplitude of the field perturbation at the beginning of the tachyonic instability, 



(66) 



As usual, this amplitude characterizes the contribution of a logarithmic interval of momenta into (<ys 2 (x)). 

2. M 2 /-y <P 2 < (e 2 - M 2 )/a-m 2 

This interval of momenta is dangerous, as the perturbations undergo the tachyonic amplification. Since the modes 
are still inside the Hubble radius, one can neglect the gravitational effects on the "frequency" co when estimating the 
growth of field perturbations. By the time the tachyonic amplification ends up, the amplitude of field perturbations 
is estimated as 



ip ~ ipi exp 



|wi|di 



Pi exp 



pi 



P3 Mdp 
H p 



(67) 



Recall that p\ ~ e/y/a and pz = M/y^y. The largest value of \uj\\ is approximately given by ([63]). Substituting this 
value into Eq. (|67|) and recalling that the background changes slowly {H ~ const), one finds that the amplitude of 
the field perturbation after exit from the tachyonic regime is of order 



exp 



a(7 + a) H 



log 



2_L 

a M 



(68) 



where we used Eq. (|66|) . 

Recall now that H is of the same order as M/y/3a during the phantom phase. Hence the large ratio e/M leads to 
a strong amplification of field perturbations. From (fT8|) . the homogeneous field <f> at the phantom and de Sitter phase 

is estimated as <j) : 
to the constraint 



^ g " pl . The requirement that the perturbation ip is smaller than the background field <f> leads 



exp 



- / 7 e i og ( 1JL 

2\Ja(j + a)H \ \ a M 



< Vs. 



-Hm 



air- 



pl 



(69) 



As an example, in the case m = H = 10 42 GeV, M = \/3aHo, a — 1 and 7 = 1/2, we obtain the constraint 
e/M < 70. As long as a and 7 are of order one, the ratio e/M should not be too much larger than unity. 
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3. < p 2 < m 2 M 2 /e 2 



After the Hubble radius crossing, the effect of the cosmic expansion can no longer be neglected when estimating 
the "frequencies" of the field perturbations. Since there are no tachyonic instabilities for m 2 M 2 /e 2 < p 2 < M 2 /^, 
we consider the evolution of perturbations in the region < p 2 < m 2 M 2 /e 2 . In Ref. [H the equations for the field 
perturbations were derived in the slow-rolling background under the condition p 2 <C M 2 ,m 2 . The equation for the 
perturbation \ = o-f is approximately given by 

d 2 /. Id 2 a 2 m 2 M 2 \ , . 

Here we neglected the contribution of metric perturbations on the r.h.s. of this equation. Note that metric pertur- 
bations works as a back reaction effect after the field perturbation is sufficiently amplified. The growth rate of the 
perturbation \ is mainly determined by the terms in the parenthesis of Eq. (|70[) rather than the backreaction of metric 
perturbations. 

The last term corresponds to the tachyonic mass term, which already appeared in Minkowski spacetime [see Eq. (|57[) ]. 
Since the term (d 2 a/dr] 2 )/a is of order a 2 H 2 , one can estimate the ratio of the tachyonic mass relative to this 
gravitational term: 

_ a 2 m 2 M 2 /e 2 ^ m 2 M 2 

= (d 2 a/d?7 2 )/a 'I 2 !! 2 '' (1) 

If we use the de Sitter value H = Mj \/3a, this ratio is estimated as S = 3am 2 /e 2 <C 1. Hence the gravitational term 
(d 2 a/dr; 2 )/a dominates over the tachyonic mass. 

In the de Sitter background with a = — 1/ (Hrj) the approximate solutions to Eq. (fTTJj) can be obtained by setting 
\ = — (C / H)r)~ 1+S . One finds that 8 — —m 2 M 2 /(3H 2 e 2 ) for the growing solution, thereby giving 

if = Crf%5& K a <5/3 _ ( 72 ) 

In Ref. [HI it was shown that the physical temporal component of the vector field perturbations evolves as bo /a oc a s ^ 3 , 
whereas the physical spatial component of the vector field decreases as Bi/a cx a~ 1+s / 3 . The growth rate of tp and 
bo /a is small due to the condition i « 1. So, the second tachyonic instability is harmless for the past and present 
cosmological evolution. However, we notice that since the de Sitter solution is a late-time attractor, the perturbations 
tp and bo/ a become larger than the homogeneous background fields in the distant future. At this stage we expect that 
the contribution of metric perturbations can be also important. 



V. CONCLUSIONS 



In this paper we have studied the dynamics of dark energy in a Lorentz-violating model with the action given in 
([1]). The model involves a vector field and a scalar field <& with mass terms M and m, respectively. The presence 
of the one-derivative term ed^&B 11 leads to an interesting dynamics at the IR scales larger than e _1 . The phantom 
equation of state can be realized without having ghosts, tachyons or superluminal modes in the UV region. 

We have taken into account the contributions of radiation and non-relativistic matter and studied the cosmological 
evolution of the system. Interestingly, there exists a de Sitter attractor solution that can be used for the late-time 
acceleration. The phantom regime is not an attractor, but we have found that in a range of parameters, the phantom 
stage occurs during the transition from the matter epoch to the final de Sitter attractor. As is seen, e.g., in Fig. [1] 
the equation of state parameter w of dark energy crosses the cosmological constant boundary towards the phantom 
region. We clarified the conditions under which the w = — 1 crossing is realized together with the existence of the 
stable de Sitter solution. 

In the model studied in this paper, the effective Newton's constant is time-dependent. We have found, however, 
that this dependence is typically mild, though for interesting values of parameters it is close to the experimental 
bounds. 

We have also considered the field perturbations in Minkowski spacetime and obtained the momentum scales of 
instabilities present in the IR region (p < e). We have found that either tachyons or ghosts appear for the spatial 
momenta p smaller than */ (e 2 — M 2 )/a — m 2 , while in the UV region there are no unhealthy modes. In the cosmo- 
logical context the presence of tachyons at the IR scales leads to the amplification of large-scale field perturbations 
whose wavelengths are roughly comparable to the present Hubble radius. There are two tachyonic regions of spatial 
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momenta in this model: (a) one is sub-horizon and its momenta are characterized by M 2 < p 2 < (e 2 — M 2 ) /a — rn 2 ; 
(b) another is super-horizon and has < p 2 < m 2 M 2 /e 2 . In the region (a) we derived the condition under which the 
perturbations always remain smaller than the homogenous fields, see Eq. (|69|l . While the existence of the phantom 
phase requires that e > M, the condition shows that e cannot be very much larger than M. Thus the allowed 
range of e is constrained to be relatively narrow. In the tachyonic region (b) the growth of the perturbations is 
estimated as ip oc a am ^ . Since the growth rate is suppressed by the factor m 2 /e 2 , this effect is negligible in the past 
and at present, though the inhomogeneities can start to dominate over the homogenous fields in the distant future. 

There are several issues yet to be understood. The presence of the tachyonic instability on sub-horizon scales may 
lead to the variation of the gravitational potential, which can be an additional source of the late-time integrated 
Saches- Wolfe effect on the CMB power spectrum. Another property of this model is the peculiar time-dependence of 
the effective Newton's constant, which may result in interesting phenomenology 

The model studied in this paper is likely to belong to a wider class of Lorentz-violating theories exhibiting the 
phantom behaviour (see Refs. [42 ] for a number of Lorentz-violating models). It would be interesting to understand 
how generic are the features we found in our particular model — late-time de Sitter attractor, transient phantom stage, 
time-dependent Newton's constant, sub-horizon tachyons, super-horizon ghosts, etc. One more direction is to modify 
our model in such a way that it would be capable of describing inflationary epoch rather than the late-time acceleration. 
Since it is known that th e sp ectra of scalar and tensor perturbations produced during the phantom inflationary phase 
are typically blue-tilted [43|], this model may give rise to some distinct features in the CMB spectrum. 
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Appendix A 

We are going to find the effective Newton's constant that determines the strength of gravitational interactions at 
distances shorter than all scales present in our model, including e _1 , M _1 , mT 1 as well as the Hubble distance. To 
this end, we neglect the last two terms in the action ([1]), and also neglect the time dependence of the background fields 
<j> and X. We also neglect the space-time curvature of the Universe, and therefore consider our model in Minkowski 
space-time. 

Let us impose the gauge hoi = 0, where h^ v is the metric perturbation about the Minkowski background. Then the 
quadratic Lagrangian for perturbations of metric, vector and scalar fields is readily calculated, 



L 



1 



1 



Xdih 



V'OO 



b + -Xh oa 



dibj Xh 



dib + -Xh 00 



(73) 



where X 2 = B 2 is the background value. Clearly the scalar field ip decouples in our approximation, so we will not 
consider it in what follows. 

By varying the quadratic action with respect to hoo and hij, one obtains (00)- and (r/)-components of the linearized 
energy-momentum tensor for perturbations (note that T^ v = —28S '/ 'Sh^J). Specifying further to scalar perturbations 



with bi = dibh and choosing conformal Newtonian gauge, /iqo = 2$, hi 



-2'Wjj, one obtains (we keep the standard 



notation for the Newtonian potential, even though the same notation was used for the original scalar field in the main 
text) 



T ° = aI(IA$ 
T) = aXd l d 1 b L - 



Ab L ) + jxn(b + x$) , 



(74) 
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where A = didi and □ = 8q — A. The field equations for 60 and &l in the absence of sources for these fields read 

n(b + x$) = 0, 

-□6 L +X($-#) =0. (75) 

Now, the longitudinal (proportional to didj) part of the (r/)-component of the Einstein equations, in the absence of 
external anisotropic stresses, gives 

$ + * = 8TrGaXb L , (76) 

while the trace part and (OO)-component are 

# + -A($ + *) = -4 7 rGaX 2 \E' - 4irG Pext , (77) 
-A* = AnG{aX 2 £\<S> - aXA6 L ) + 4TtGp ext , (78) 

where p ox t and p e xt are energy density and pressure of an external source. 

For time-independent, pressureless source it is consistent to take all perturbations independent of time and set 
bz = 0. Then one finds, as usual, \& = — $ and obtains the following equation for the Newtonian potential, 

(1 - 47rGaX 2 )A$ = 4nGp exi ■ (79) 

Thus, the effective Newton's constant in the background field X is 

G* = G(l - AnGaX 2 )- 1 . (80) 

This means that the effective Planck mass entering the Newton's law is given by (J6j> . 

Appendix B 

A. Autonomous equations 

Let us define the following dimensionless variables which are convenient for studying the dynamical system [H, [38| : 
Ky^jX k6 Kind) [4tt X M K.Jp~r 

81 = W X2 = 7M> X3 = VM> ^VT^' X5= H> Xe =VM- (81) 

Then we obtain the following autonomous equations 

, 3a el 1 2 H' , 

x 1 = -Sxi — Xi + ——-x 2 x 5 + —-x i x 5 - xi— , (82 

n/t M Vi \ff H 

.el e m H' , 

Xn = —oxo x-ixr, — 3 — xaXc. xhXk — Xi — , (83 

2 Z M^ M M E v ' 

, m H' 

x 3 = —x 2 x 5 - x 3 — , (84) 
,1 

x i = —xi , (85) 

a 1 

x' 5 = ~x^—, (86) 
H 1 

x' 6 = -2x 6 - x 6 — , (87) 

and 

H' 31 + ^+12-^+ (3a + x 2 ) + z§/3 + 2(e/M)x 2 X4,x 5 + 4(a/ y /j)x 1 x 4 
IT ~ ~2 1 + iax 2 A ' 
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where prime denotes the derivative with respect to 

N = ln(a) . 

Equation ITT|) gives the constraint 

2 

= -j^j- = \-x\~x\-x\+x\(^ct + xf) - xl . 

Note that the above equations are invariant under the simultaneous change of the signs of <fi and X. Hence it is not 
restrictive to study the case of positive (f>. Note also that we study the case of an expanding Universe with H > 0. 



B. Fixed points 

By setting x[ = one formally finds the following six fixed points: 



(A) de Sitter (i): (x lt x 2 , x 3 , x 4 , x 5) x 6 ) = (o, 0, ^ ^Z^-^- , - ^^^^ , y/te, oj 



• (B) de Sitter (ii): (xi, X2,X3, X4, X5, x 6 ) = (0, 0, const, 0, 0, 0) , 

• (C) matter: (xi,X2,X3,X4, x^,Xe) — (0,0,0,0,0,0), 

• (D) radiation: (xi, x 2 , x 3 , X4, X5, xq) = (0, 0, 0, 0, 0, 1) , 

• (El) kinetic point (i): (xi, x 2 , X3, X4, X5, xq) = (0, 1, 0, 0, 0, 0) . 



(E2) kinetic point (ii): (x\,x 2 , X3, X4, x§,xq) = (0, — 1, 0, 0, 0, 0) . 



The fixed point (A) is precisely the de Sitter solution (|18p that we discussed in Sec. IIII Al We will comment on its 
stability shortly. 

The point (B) is also in some sense a de Sitter point. It exists even in the absence of the field X and satisfies 
the relation 3H 2 = K 2 W((f>). To reach the solution (B), the Hubble parameter needs to increase towards infinity 
[M/H — ► 0), and the field 4> needs to diverge as well. 

The point (C) corresponds to matter-dominated era satisfying fi m = 1 and w c $ = 0, whereas the point (D) describes 
radiation-dominated epoch with Q r = 1 and w c s — 1/3. 

The points (El) and (E2) are kinetic solutions satisfying f2oE = 1 and w c g — 1. These solutions are used neither 
for dark energy nor for radiation/matter dominated epochs. 

A cosmologically viable trajectory starts from the radiation point (D), connects to the matter solution (C) and 
finally approaches the de Sitter point (A). [Note that the initial data (|22p indeed correspond to xi, x 2 , X3, X4, x§ — > 
as t — > 0.] To see that this sequence of events is indeed possible, let us study the stability of the fixed points against 
perturbations. 

Let us consider linear perturbations Sxi. By perturbing Eqs. (j52"|) - (f8T>|) we obtain 
5x[ = ( -3 - — - cixi ) 5xi + ( -^--^—x 5 - c 2 xi J Sx 2 - 03X18x3 



H J \M yft 

— Xg - - C4X1 ] 5x4 + (77^12 + -7=x 4 x 5 - C5X1 ] 5x 5 - c 6 xi5x 6 , (89) 
5x 2 = - (j^-^=x 5 + cix 2 ^j 5x\ - + c 2 x 2 + 5x 2 - (j^x 5 + c 3 x 2 ^j 5x 3 

- ( 3 ]^ a; 5 + C4X2) 5x4 - \Tj~F Xl + ^~m Xa + Jj X3 + C5X2 j Sx5 ~ °e x 2Sx 6 , (90) 



5x' 3 



(771 \ f 7TI \ 

—x 5 - c 2 x 3 J 5x 2 - c 3 x 3 5x 3 - c 4 x 3 5x 4 + yr^x 2 - c 5 x 3 J 5x 5 - c 6 x 3 5x 6 , (91) 

5x' 4 = —Sxi, (92) 

5x' 5 = —CiX^Sxi — C2X5<5X2 — c 3 x$5x 3 — c 4 x§5x 4 — + c 5 x 5^j $ x 5 ~ CqX^Sxq , (93) 

5x' 6 = —c\x & 5xi — c 2 x 6 5x 2 — c 3 x 6 5x 3 — c 4 x & 5x 4 — c 5 x 6 5x 5 — ( 2 + — + c 6 x 6 J 5x 6 , (94) 
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where S(H' / H) = Y^i=i c $ x i with 



3xi + 6(a/ ^/7)x4 3a; 2 + 3(e/M)x 4 a;5 3x 3 

C l = 2 X4 , C 2 = — — g ! C 3 = 771 2 ' 

1 + 6ax\ 1 + 6ax\ 1 + 3aa;| 

3x4(3a + x\) + 3(e/M)x2Xs + 6(a/y/j)xx 6ax± H' 
1 + iaxl 1 + 2>ax\ H 

3x^X5 + 3(e/M)x2X4 x 6 

1 ~t~ 3q^4 



c 4 



^ ""*4"<o i "V 1 -/ y^z^^i ^ (95) 

1 + 'icxx 1 , 1 + 3ax? 



The stability of fixedpoints can be analyzed by considering eigenvalues of the 6x6 matrix M. for perturbations 
along the lines of Ref. [3 [H[ . 

The stability of the de Sitter point (A) is important for having the late-time accelerated epoch. This depends upon 
the two ratios e/m and M/m once the parameters a and 7 are fixed. When a — 1 and 7 = 1/2, for example, the 
parameter range of e/m is determined by the ratio M/m. We find that the point (A) is a stable attractor if the 
following conditions hold: 

• (i) When M/m = 0.1, e/m > 0.817, 

• (ii) When M/m = 1, e/m > 1.35, 

• (hi) When M/m = 10, e/m > 3.52, 

When f»ra, the stability of the point (A) is ensured automatically unless the ratio M/m is too much larger than 
unity. In view of (f20|) the case e 3> m is of particular interest. 
For another de Sitter point (B), the eigenvalues are 



-3, -3,-|±i^9-— ,0,-1/2, (96) 

This means that this point is marginally stable. The zero eigenvalue comes from the perturbation equation for 5x$. 
If H continues to increase toward the solution (B), this eigenvalue actually obtains a small negative value, as can be 
seen from Eq. (|93| . Thus in such a case the point (B) is stable. However, we know that the phantom phase is realized 
only for finite field values bounded by 4>a, see Eq. (|34|) . Hence it is not possible that the actual solutions approach 
the point (B) with infinite H and <fi. 

The matter point (C) has the eigenvalues 



4 / 48a 

3/2, -3/2, --± ^9- -J^, 0,-1/2, (97) 

which shows that the matter era corresponds to a saddle point with one positive eigenvalue. Hence the solutions 
eventually repel away from this fixed point even if they temporarily approach it. 
The radiation point (D) has the eigenvalues 



-ihsV 1 - 1 ^- 1 ' ' 2 ' 1 ' (98) 

One finds that the radiation epoch corresponds to a saddle point with two positive eigenvalues. 
The kinetic points (El) and (E2) have the eigenvalues 

3,3,1,0,-^5,-^5, (99) 

which shows that they are unstable. 

The above stability analysis shows that the sequence of radiation, matter and de Sitter epochs can indeed be 
realized. 
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